A chiral quark model study of Z + (4430) in the molecular picture 
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We investigated the bound state problem of the S wave charged D\D* (D[D*) system in a chiral 
quark model by solving the resonating group method equation. Our preliminary study does not 
favor the molecular assumption of Z + (4430). On the contrary, if Z + (4430) is really a molecule, 
its partner with opposite G-parity should also exist and probably may be found in the Tv + r] c (2S), 
J/ipn + iT , or tp'n + n° channel. For the bottom systems, we found the existence of both I°=l + and 
I a =l~ BiB* (B[B*) molecules is possible. 

PACS numbers: 12.39.-x, 12.40.Yx, 13.75.Lb 



I. INTRODUCTION 

In recent years, a series of heavy quark hadrons with 
unexpected properties were observed one by one, from 
D s j(2317) [J & H, X(3872) LJEJEH], Y(4260) 
@, II Gil, V(3940) EL and Y(3940) ©El, to the newly 
observed y(4140) [14]]. These near threshold mesons 
stimulated the interpretations beyond the quark model. 
For the interesting hidden charm XYZ states, the inter- 
pretations include tetraquark or molecular states, hybrid 
charmonia, cusps, and threshold effects. However, it is 
not excluded that these XYZ states are still mesons dom- 
inated by cc components. On the contrary, the observa- 
tion of charged charmonium-like states is a surprising 
issue because such states contain at least four quarks. 

The Belle Collaboration announced a distinct peak 
i? + (4430) in the 7r + ^' invariant mass distribution in 
the decay B — > K^tp' in Ref. [HI]. The mass and 
width are M = 4433 ± A(stat) ± 2(syst) MeV and 
r = 45±\l(stat)±H(syst) MeV, respectively. The mini- 
mum quark content is ccud. Very recently, a little heav- 
ier and broader Z + (4430) is obtained in Belle's reanalysis 
based on the same data sample [l~6j |. However, the experi- 
mental data from the BaBar collaboration do not provide 
significant evidence for the existence of Z + (4430) [171 ] . 

In addition to Z + (4430), Belle Collaboration recently 
observed two more charged resonance-like structures in 
the 7r + y e i invariant mass distribution in B — > Kn^ip 1 
decays [18[. The mass and width for the first structure 
are 



Mi = 4051 ± 14(stat)±|?(syst) MeV 
Tj = 82i^(stat)±^(syst)MeV 

while the values for the second one are 



M 2 = 4248±^(stat)±if (syst)MeV 
T 2 = 177lg(stat)t^ 6 (syst)MeV. 



(1) 

(2) 

(3) 
(4) 



Whether the existence of these two structures is sup- 
ported awaits experimental confirmation from other col- 
laborations. 

Since the announcement of Z + (4430), lots of discus- 
sions in various pictures have appeared, which include a 
tetraquark state [2(1 . a resonance or a molecule in the 
D X D* [DyD*) channel [U 1| H , a baryonium 
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26|, a threshold cusp [2jJ and a radial excited cs 
28[. In Ref. the bottom analogs of Z+(4430) 
were studied. Besides the spectroscopy, there were dis- 
cussions about its production [3(1 [U, H3, [H, l3~j and 
decay [35|]. The tti/j' scattering is studied in Ref. [3a ]. 

The molecular picture is widely used because Z + (4430) 
is close to the threshold of D[D* or D\D*. In Ref. [H, 
the calculation at hadron level indicates that it is possi- 
ble to get a bound state in the D[D* or D±D* system 
with appropriate parameters and to interpret Z + (4430) 
as a molecule. The QCD sum rule study (24| and a quark 
model calculation [25|] also favor the D\D* molecule in- 
terpretation. However, a recent calculation on the lattice 
indicates such an interpretation is probably problematic 

ma 

To help to understand this charged state further, we 
present our preliminary study from a chiral quark model 
(yQM) [3l| and an extended chiral quark model (E%QM) 
[391 ] calculation in this article. The former model includes 
a and tt exchange interactions between light quarks. The 
later model is an extended version of the former one by 
including the vector meson exchanges. We investigate 
the bound state problem of the S-wave D[D* or D\D* 
system by solving the resonating group method (RGM) 
equation [4(J. In previous studies, this approach has 
successfully reproduced the energies of the light quark 
baryon states, the binding energy of the deuteron and 
the NN scattering phase shifts. When using it to study 
the system of a light meson and a light baryon [4l|, the 
resulting phase shifts are also in agreement with the ex- 
perimental data. With this model, we have preliminarily 
studied the bound state problem of two S-wave heavy 
mesons in Ref. (42[ and [H|. The results are rou ghly 
consistent with similar studies at hadron level 0, |45| . 
We here intend to explore whether or not the model can 
be used to the case of orbitally excited mesons. 
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For the system studied, the orbitally excited heavy 
mesons are Di and D[. These two J p =l + mesons are 
mixed from 3 Pi and 1 P\ states 



\D X ) = cos0| x Pi) +sin6»| 3 Pi) 
\D[) = -s33x6\ x Pi) +cos6»| 3 Pi) 



(5) 



The mixing angle 9 = —54.7° or 35.3° may be de- 
duced with the mass of the heavy quark going into in- 
finity. In this article, we adopt the widely adopted value 
9 = -54.7° [HHI. 

In the molecular picture, the flavor wave function of 
Z+(4430) reads 



1 



Z+(4430) = — -=(D°D* 
v2 



D *o D +\ 



1 



Z+(4430) = -^={D'?D* + + D*"D'f) 



s*o n'+\ 



(6) 



(7) 



The quantum numbers are I G J p —l + (0, 1, 2)~. As a pre- 
liminary study, we consider only interactions involving 
color-singlet mesons. The one-meson exchange potentials 
between the heavy mesons are induced by the meson ex- 
changes between light quarks. 

This paper is organized as follows. After the introduc- 
tion, we present our model in Sec. |TTJ Then in Sec IIII1 
we present our results and discussions. 



II. THE CHIRAL QUARK MODEL 

The Hamiltonian for the heavy quark meson-antimeson 
system in the chiral quark model has the form [H, [39[ 



V 



OGE 



V 



conf 



M 



(8) 



where Tj is the kinetic term of the ith quark or antiquark 
and Tg is the kinetic energy operator of the center of 
mass motion. M is the exchanged meson between light 
quarks. In the chiral quark model, one of the sources 
for the constitute quark mass is the coupling with chiral 
fields which come from the spontaneous vacuum break- 
ing. Because the breaking has small effects on the gen- 
eration of the constitute mass of the heavy quarks, the 
coupling of the a meson and the heavy quarks should be 
weak. As a result, the possible flavor-singlet meson ex- 
change interactions between heavy quarks and between a 
heavy quark and a light quark have small contributions 
and so we ignore them. 

The potential induced by the one-gluon-exchange 



(OGE) interaction reads 



V, 



OGE 



1 7T 



1Q 



r 2 
4 1 / 



1 1 



3 m q mQ 



' v OGEi 



(9) 



V, 



1 



OGE 



j for quarks and - ■ 



A for antiquarks 



where F; 

and m q {tuq) is the light (heavy) quark mass. The linear 
confinement potential inside the color-singlet meson is 



V qQ 



'q- r Q \ a qQ T f ifo) 



There are similar expressions for V^q }e and V^q 1 . Be- 
cause we preliminarily ignore the possible hidden color 
contributions, we do not need V^q E and V qq GE ■ 

For the meson exchange potentials between two light 
quarks, we have [38l. |39| 



-C{g chl m<j, A)X 1 {m„, A, n 



V w "{ r ij) = C(g ch ,m^ a ,A)— — ^ — X 2 (m„ a , A, r y ) 



12m qi m qj 



x[er(i) ■ tr(j)][T a (i)r a (j)], 
(a = 1,2,3) 



(11) 



(12) 



rrr 



C(g chv ,m Pa ,A) ^ X 1 (m Pa ,A : r ij ) + p — — 



X 1 + 



fchv m qi + m qj ^ / fchv \ 2 m H W 9j 
9chv Mn 9chv M N 

xX 2 (m Pa , A, nj 

)kM ■<?)]} [Ta(i)T a (j)}, 



(13) 



v uu ( r v ) = C(9chv , m u , A) <^ Xi (m u , A, ry ) + 



6m: 



g c hv Mn g c hv M N 



xX 



<.(m u ,A,rij)[tr(i) ■ tr(j)] j 



(14) 



where 



C(g ch ,m,A) = 



9ch A2 ™ 



A. 



Xi(m,A,r) = Y(mr) Y(Ar) 



(15) 
(16) 



X 2 {m,A,r) = Y(mr) - I — ) Y(Ar), (17) 



Y(x) = — . 



(18) 



We do not present the tensor term and the spin-orbital 
term in the meson exchange potentials since we consider 
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XQM E X QM 

Set 1 Set 2 Set 3 

b u (fm) 05 0.45 0.45 
m u (MeV) 313 313 313 
m a (MeV) 595 535 547 
g chv 2.351 1.972 

fchv/gchv 2/3 

TABLE I: Three sets of model parameters. Other meson 
masses are: = 138 MeV, m p — 775.8 MeV, and = 
782.6 MeV. 



only S-wave meson-meson interactions. Here we use the 
same cutoff A in describing various meson interactions. 
Its value is around the scale of chiral symmetry breaking 
(-1 GeV). 

The interaction between a quark and an antiquark is 
related to that between two quarks through the rela- 
tion VM = GuVqq where Gm is the G-parity of the 
exchanged meson. 

By calculating the RGM matrix elements and solving 
the RGM equation for the bound state problem, one gets 
the energy of the system and the relative motion wave 
function, from which one deduces the binding energy 
E n = Mq q + MQq — M sys . If E is positive, the sys- 
tem is bound. 

For the model parameters, we take the values de- 
termined in the previous investigations [H, H^]. The 
harmonic-oscillator width parameter b u =0.5 fm for xQM 
and b u =0A5 for E%QM. The up (down) quark mass 
™«(d)=313 MeV. The coupling constant g c h— 2.621 is 
derived from the measured NNt: coupling constant 
Snn-k/^ 71 = 13.67. The masses of 7r, p, and io are taken 
to be the experimental values, whereas a meson mass 
is adjusted to fit the binding energy of the deuteron. 
In ExQM, we use two sets of values. We present the 
above parameters in Table HI The parameters in the 
OGE and the confinement potentials can be derived from 
the masses of the ground state baryons and the heavy 
mesons. In fact, their values do not give effects to the 
binding energy of the meson-antimeson system when we 
ignore the hidden color contributions [42] . So we do not 
present their values here. For the charm quark masses, 
we take m c =1430 MeV H3 and 1870 MeV [H] to see the 
heavy quark mass dependence of the binding energy. For 
bottom quark, we use m & =4720 MeV g| and 5259 MeV 
48]. In our calculation, we take two values for the cutoff 
A=1000 MeV and A=1500 MeV. 



III. RESULTS AND DISCUSSIONS 

We first study the S-wave D\D* system. We illus- 
trate the diagonal meson exchange matrix elements of the 
Hamiltonian in the generator coordinate method (GCM) 
calculation for different angular momenta in Fig. [T] 
with the parameters Set 3 in Table HI m c =1870 MeV 
and A=1500 MeV. One finds the dominate contributions 
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(c) 

FIG. 1: The meson exchange GCM matrix elements for J=0 
(a), J=l (b), and J=2 (c) D\D* system. The used parameters 
are 6„=0.45 fm, ma=547 MeV, g chv =i-972, fchv/gchv^/3, 
m c =1870 MeV, and A=1500 MeV. 
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come from the a and ir meson exchange interactions. 
The J=0 system is more attractive than J—l, 2 systems. 
However, we do not find a binding solution in this system 
with the parameters presented in the former section. For 
the S-wave D[D* system, the meson exchange GCM ma- 
trix elements are illustrated in Fig. [2] The J— 2 system 
is the most attractive one. But the system is also un- 
bound. So our preliminary calculation does not support 
the interpretation that Z + (4430) is an S-wave D\D* or 
D[D* bound state. _ 

The D X D* or D[D* can also form a G=- system with 
the flavor wave function 



z' = 



D* U D+) 



(19) 



(a) 



Z' 



V2 



'o n*+ 



(D»D 



*o n '+\ 



D* U D 



(20) 



We found these systems are also unbound with our pa- 
rameters. The attractive force in this case is a little 
stronger than that in the G=+ case. One observes this 
feature by comparing the GCM matrix elements for J=Q 
case in Fig. [3] with those in diagram (a) of Fig. Q] 

The bottom analogs have better chances to form 
molecular states because the kinetic term in the Hamilto- 
nian has relatively small contributions. We do get bind- 
ing solutions with the parameters in the former section. 
Table [II] gives the binding energy and the root-mean- 
square (RMS) radius for the B\B* and B[B* systems. 
The system is unbound for J=2 B X B* and J=0, 1 B[B* . 
In that table, we present both the results for the G=+ 
case and those for the G=— case. A little deeper bound 
states appear in the later case. 

Our study with the chiral quark model approach does 
not support the existence of an S wave molecule in the 
D\D* and D[D* systems. This result is inconsistent 
with our similar study at hadron level [23| . However, in 
the case of DD* system, we got consistent conclusions 
with these two approaches [42|, |44| . A possible reason 
for the present inconsistency is due to the different ap- 
proximations in getting the potentials. As a first step 
to derive the potential, one writes out the quark-quark 
(or meson-meson) scattering matrix in momentum space. 
The denominator of the propagator for a meson reads 
p 2 — m 2 + ie — Pq — p 2 — m 2 + ie where p (p) is the 
four(three)-momentum and m is the meson mass. The 
approximation p 2 — > — p 2 , i.e. po ~ 0, is adopted in the 
chiral quark model approach, whereas the possible large 
Po is considered for the hadron level calculation [23|. In 
the later approach, the principal integration is always 
assumed if po is larger than the meson mass m when 
we get the coordinate-space potentials. In the present 
case, the large p<j is around 3m T while the large p<j is 
about m w +7 MeV in the DD* case. Probably it is this 
Po around Sm^ leads to inconsistent conclusions for the 
studies using these two approaches. We reanalyzed the 
binding energies of the D\D* (D[D*) system at hadron 



S (fir 



(b) 



(c) 

FIG. 2: The meson exchange GCM matrix elements for J=0 
(a), J=l (b), and J=2 (c) D[D* system. The used parameters 
are the same as those for D\D* . 
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TABLE II: The binding energy (RMS radius) for the bottom analog of DiD* 
means the system is unbound. The system is unbound for J=2 BiB* and J=0 



and D[D* system in unit of MeV (fm). 
1 B^B*. 



A x 



G-parity m b (MeV) A (MeV) 
G=+ 4720 1000 
1500 

5259 1000 
1500 

G=- 4720 1000 
1500 

5259 1000 
1500 



B 1 B* (J=0) 
Set 1 Set 2 Set 3 

1.0(1.4) 3.7(1.2) 3.4(1.2) 

1.6(1.4) 4.8(1.2) 4.4(1.2) 

1.8(1.3) 4.8(1.1) 4.5(1.1) 

2.4(1.3) 5.9(1.1) 5.5(1.1) 

5.2(1.1) 9.9(1.0) 9.4(1.0) 

6.9(1.1) 13.0(0.9) 12.4(0.9) 

6.5(1.1) 11.8(0.9) 11.3(0.9) 

8.4(1.0 15.1(0.9) 14.5(0.9) 



SiB* (J=l) 

Set 1 Set 2 Set 3 

x 1.2(1.3) 0.9(1.3) 

x 1.9(1.3) 1.6(1.3) 

x 2.0(1.2) 1.7(1.2) 

0.4(1.4) 2.9(1.2) 2.5(1.2) 

1.3(1.3) 4.3(1.1) 3.9(1.1) 

2.2(1.3) 6.0(1.1) 5.5(1.1) 

2.3(1.2) 5.6(1.1) 5.2(1.1) 

3.3(1.2) 7.5(1.0 7.0(1.0 



B[B* (J=2) 

Set 1 Set 2 Set 3 

x 0.6(1.3) 0.3(1.3) 

x 1.3(1.3) 1.0(1.3) 

x 1.4(1.2) 1.1(1.2) 

x 2.3(1.2) 1.9(1.2) 

x 1.3(1.3) 0.9(1.3) 

x 2.2(1.2) 1.9(1.2) 

x 2.2(1.2) 1.9(1.2) 

0.6(1.3) 3.3(1.1) 2.9(1.1) 



50.0 A 1 1 1 

40.0- '~'--v % | 1 

3o.o- \ * 

20.0- 
p 10.0- 

>~10.0- '/^^^ 

-30.0- J' y' 

-30.0- S?y'' 
-40 - ' 

-50.0 \ 1 1 1 

0.0 0.5 1.0 1.5 2.0 

S (fm) 

FIG. 3: The meson exchange GCM matrix elements for J=0 
D\D* system with G—-. The used parameters are the same 
as in Fig. [T] 



level with the approximation po ~ 0. As expected, we 
did not find a binding solution, which indicates the im- 
portant role of po- However, we need the experiments to 
judge which approximation is correct. The comparison 
of model predications with experimental measurements 
may finally answer the puzzle. 

Although our model calculation does not support the 
assumption that Z' + (4430) is a molecule, such an inter- 
pretation is still possible. To get a more conclusive result 
in a future investigation, the following effects may be in- 
cluded. First, the hidden-color configuration and a larger 
model space may have contributions and can be consid- 
ered. Secondly, the coupling with D wave interaction is 
probably not negligible and may be studied. Thirdly, the 
different approximation in deriving the coordinate space 



potential may be investigated in detail. In addition, our 
model neglects the contribution from a exchange interac- 
tion between two heavy quarks or between a heavy quark 
and a light quark. Although the coupling constant gQQ a 
is expected to be small, the value may have big effects 
because no mass factor in the a potential can suppress 
the contribution. This is also an open question one may 
discuss. 

According to the GCM matrix elements, which roughly 
reflect the force between the two mesons, if Z + (4430) can 
be identified as a D\D* or D[D* molecular state, a G=- 
state around 4430 MeV should also exist. One expects 
that such a state may be searched for in the ir + j] c (2S), 
J/ifjir + ir°, or tp'ir + ir° channel. 

In short summary, we have studied the bound state 
problem of the S wave D\D* (D^D*) system in a chiral 
quark model. Our preliminary calculation does not favor 
the assumption that Z + (4430) is an S wave molecule. On 
the contrary, once Z + (4430) (G=+) may be identified as 
a D 1 D* (D[D*) molecule, its partner with G = - should 
also exist. When we move on to the bottom analogs, the 
existence of the charged B\B* (B[B*) molecules with 
G=+ and G=— are both possible. Such states can prob- 
ably be found in the nT{2S), 7^(25), T(15)7r+7r°, and 
T(lS)ir + ir° channels in future measurements. 
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